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> ! Abstract 
O 

' We present in detail the calculation of the virtual 0(a s ) corrections to the 

inclusive semi-leptonic rare decay b — > si + i~ . We also include those 0(a s ) 
bremsstrahlung contributions which cancel the infrared and mass singularities 
showing up in the virtual corrections. In order to avoid large resonant con- 
Q_i! tributions, we restrict the invariant mass squared s of the lepton pair to the 

p ' range 0.05 < s/mf < 0.25. The analytic results are represented as expansions 

in the small parameters s = s/mf, z = rn^/mf and s/(4m^). The new contri- 
butions drastically reduce the renormalization scale dependence of the decay 
spectrum. For the corresponding branching ratio (restricted to the above s- 
range) the renormalization scale uncertainty gets reduced from ~ ±13% to 
S 1 ~ ±6.5%. 



*Work partially supported by Schweizerischer Nationalfonds and SCOPES program 



I. INTRODUCTION 



Rare -B-decays are an extremely helpful tool for examining the standard model (SM) and 
searching for new physics. Within the SM, they provide checks on the one-loop structure 
of the theory and allow one to retrieve information on the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements V ts and V t d, which cannot be measured directly 

The first measurement of the exclusive rare decay B — > K*^ was obtained in 1992 by 
the CLEO collaboration Somewhat later, also the inclusive transition B — > X s ^ was 
observed by the same collaboration Although challenging for the experimentalists, the 
inclusive decays are clean from the theoretical point of view, as they are well approximated 
by the underlying partonic transitions, up to small and calculable power corrections which 
start at 0(A 2 QCD /ml) |,|. 

The measured photon energy spectrum [0 and the branching ratio for the decay B — > X s 7 
0J^J^] are in good agreement with the next-to-leading logarithmic (NLL) standard model 
predictions (see e.g. P|-fH{|). Consequently, the decay B — > X s j places stringent constraints 



on the extensions of the SM, such as two-Higgs doublet models [|T~Q| , p~5| , p~6|1 , supersymmetric 



models [p^ -|22|], etc 



B — > X s £ + £~ is another interesting rare decay mode which has been extensively consid- 
ered in the literature in the framework of the SM and its extensions (see e.g |^-|28|). This 
decay has not been observed so far, but it is expected to be measured at the operating B- 
factories after a few years of data taking (for upper limits on its branching ratio we refer to 
1 29| , |30[] ) . The measurement of various kinematical distributions of the decay B — > X s £ + £~, 



combined with improved data on B — > X s 7, will tighten the constraints on the extensions 
of the SM or perhaps even reveal some deviations. 

The main problem of the theoretical description of B — > X s £ + £~ is due to the long- 
distance contributions from cc resonant states. When the invariant mass ^/s of the lepton 
pair is close to the mass of a resonance, only model dependent predictions for such long 
distance contributions are available today. It is therefore unclear whether the theoretical 



uncertainty can be reduced to less than ±20% when integrating over these domains 31]. 

However, restricting -y/i to a region below the resonances, the long distance effects are 
under control. The corrections to the pure perturbative picture can be analyzed within the 
heavy quark effective theory (HQET). In particular, all available studies indicate that for 
the region 0.05 < s = s/ml < 0.25 the non-perturbative effects are below 10% |32]-|37|. 
Consequently, the differential decay rate for B — > X s £ + £~ can be precisely predicted in 
this region using renormalization group improved perturbation theory. It was pointed out 
in the literature that the differential decay rate and the forward-backward asymmetry are 
particularly sensitive to new physics in this kinematical window |38|-fi0 . 



Calculations of the next-to-leading logarithmic (NLL) corrections to the process B — > 
X s £ + £~ have been performed in refs. |2l] and [^j . It turned out that the NLL result suffers 



from a relatively large (±16%) dependence on the matching scale nw To reduce it, next- 
to-next-to leading (NNLL) corrections to the Wilson coefficients were recently calculated by 
Bobeth et al. ETJ. This required a two-loop matching calculation of the effective theory to 
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the full SM theory, followed by a renormalization group evolution of the Wilson coefficients, 



using up to three- loop anomalous dimensions |41JTll. Including these NNLL corrections to 



the Wilson coefficients, the matching scale dependence is indeed removed to a large extent. 

As pointed out in ref. [M, this partially NNLL result suffers from a relatively large 
(~ ±13%) renormalization scale (//&) dependence (/if, ~ 0{m b )) which, interestingly enough, 
is even larger than that of the pure NLL result. Recently we showed in a letter f42| that 
the NNLL corrections to the matrix elements of the effective Hamiltonian drastically reduce 
the renormalization scale dependence. The aim of the current paper is to present a detailed 
description of the rather involved calculations and to extend the phenomenological part. 
We will discuss in particular the methods which allowed us to tackle with the most involved 
part, viz. the calculation of the 0(a s ) two-loop virtual corrections to the matrix elements 
of the operators 0\ and 0<i- We also comment on the 0(a s ) one-loop corrections to O-j-Oiq. 
Furthermore, we include those bremsstrahlung contributions which are needed to cancel 



infrared and collinear singularities in the virtual corrections. As shown already in [E], the 



new contributions reduce the renormalization scale dependence from ~ ±13% to ~ ±6.5%. 

The paper is organized as follows: In chapter [H] we review the theoretical framework. Our 
results for the virtual 0(a s ) corrections to the matrix elements of the operators 0\ and 2 
are presented in chapter III, whereas the corresponding corrections to the matrix elements 
of O7, 0$, Og and 0\q are given in chapter IV. Chapter V is devoted to the bremsstrahlung 
corrections. The combined corrections (virtual and bremsstrahlung) to b — > s£ + £~ are 
discussed in chapter VI. Finally, in chapter VII, we analyze the invariant mass distribution 
of the lepton pair in the range 0.05 < s < 0.25. 



II. EFFECTIVE HAMILTONIAN 



The appropriate framework for studying QCD corrections to rare -B-decays in a system- 
atic way is the effective Hamiltonian technique. For the specific decay channels b —>■ s£ + £~ 
{£ = /i, e), the effective Hamiltonian is derived by integrating out the heavy degrees of 
freedom. In the context of the standard model, these are the t-quark, the H^-boson and the 
Z°-boson. Due to the unitarity of the CKM matrix, the CKM structure factorizes when 
neglecting the combination V* s V u b. The effective Hamiltonian then reads 

4GV A 

«eff = --7f V tsV tb £ CiQl)O t QM) . (1) 
* 1=1 
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Following ref. fLLj , we choose the operator basis as follows: 



Ox = 


/ — rnn \ f — / / rnn 7 \ 

(s L -f^T a c L )(c L ^T a b L ) , 


o 2 = 


/ — \ f — 11I \ 

{sLl^c L ){c L ^b L ) , 


o 3 = 


(sz7A) E g (97 M ?) , 


4 = 


(s L l»T a b L ) EgiQl^q) , 


o 5 = 


{sLl^ulpbL) EgCgVW?) , 


o 6 = 


{sLl^l P T a b L ) E q (<lYYYT a q) , 


o 7 = 


^m b (s L( r^b R )F^ , 


Os = 


j-m b (s L a^T% R )G% , 


o 9 = 


g(^7A)E/(^)> 


O10 = 


§ fe7A) £#7^750, 



(2) 



where the subscripts L and R refer to left- and right- handed components of the fermion 
fields. 

The factors 1/gl in the definition of the operators O7, O9 and O10, as well as the factor 
l/g s present in 0§ have been chosen by Misiak |24j] in order to simplify the organization of 



the calculation: With these definitions, the one-loop anomalous dimensions (needed for a 
leading logarithmic (LL) calculation) of the operators are all proportional to g 2 , while 
two-loop anomalous dimensions (needed for a next-to-leading logarithmic (NLL) calculation) 
are proportional to g*, etc.. 

After this important remark we now outline the principal steps which lead to a LL, NLL, 
NNLL prediction for the decay amplitude for b — > s£ + £~: 

1. A matching calculation between the full SM theory and the effective theory has to be 
performed in order to determine the Wilson coefficients Cj at the high scale fiw ~ 
mw,m t . At this scale, the coefficients can be worked out in fixed order perturbation 
theory, i.e. they can be expanded in g 2 s : 

cm = cf W) + Sr^ M + jS^ c ^ M + ° {g " s) ■ (3) 

At LL order, only Cf' are needed, at NLL order also C\ l \ etc.. While the coefficient 
Cj 2 \ which is needed for a NNLL analysis, is known for quite some time §, Cf } and 
C{q have been calculated only recently |y] (see also |j43"|). 



2. The renormalization group equation (RGE) has to be solved in order to get the Wilson 
coefficients at the low scale fi b ~ mj. For this RGE step the anomalous dimension 
matrix to the relevant order in g s is required, as described above. After these two 
steps one can decompose the Wilson coefficients Ci(fib) into a LL, NLL and NNLL 
part according to 

CM = cfW + f^fV) + §^ c ? ] ^) + °^s) • ( 4 ) 

3. In order to get the decay amplitude, the matrix elements | (/xj,) 1 6) have to be 
calculated. At LL precision, only the operator Og contributes, as this operator is the 
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only one which at the same time has a Wilson coefficient starting at lowest order and 
an explicit l/g% factor in the definition. Hence, at NLL precision, QCD corrections 
(virtual and bremsstrahlung) to the matrix element of Og are needed. They have 



been calculated a few years ago p4|j2q| . At NLL precision, also the other operators 
start contributing, viz. Oj(fib) and Oio(a*&) contribute at tree-level and the four-quark 
operators 0\,...,Oq at one-loop level. Accordingly, QCD corrections to the latter 
matrix elements are needed for a NNLL prediction of the decay amplitude. 

The formally leading term ~ (1 / g^)Cg°\/j,b) to the amplitude for b — > s£ + £~ is smaller 
than the NLL term ~ (l/g 2 s ) [g 2 s / (16n 2 )] c£\fi b ) ||. We adapt our systematics to the 
numerical situation and treat the sum of these two terms as a NLL contribution. This is, 
admittedly some abuse of language, because the decay amplitude then starts out with a 
term which is called NLL. 

As pointed out in step 3), 0(a s ) QCD corrections to the matrix elements (s£ + £~\Oi(fib)\b) 
have to be calculated in order to obtain the NNLL prediction for the decay amplitude. In 
the present paper we systematically evaluate virtual corrections of order a s to the matrix 
elements of Ox, 2 , 7 , O s , 9 and Oi . As the Wilson coefficients of the gluonic penguin 
operators 0^,...,Oq are much smaller than those of Ox and 2 , we neglect QCD correc- 
tions to their matrix elements. As discussed in more detail later, we also include those 
bremsstrahlung diagrams which are needed to cancel infrared and collinear singularities 
from the virtual contributions. The complete bremsstrahlung corrections, i.e. all the finite 
parts, will be given elsewhere jPJ. We anticipate that the QCD corrections calculated in 



the present paper substantially reduce the scale dependence of the NLL result. 



III. VIRTUAL 0{as) CORRECTIONS TO THE CURRENT-CURRENT 

OPERATORS Ox AND 2 

In this chapter we present a detailed calculation of the virtual 0(a s ) corrections to the 
matrix elements of the current- current operators Ox and 2 . Using the naive dimensional 
regularization scheme (NDR) in d = 4 — 2e dimensions, both, ultraviolet and infrared singu- 
larities show up as l/e n poles (n = 1,2). The ultraviolet singularities cancel after including 
the counterterms. Collinear singularities are regularized by retaining a finite strange quark 
mass m s . They are cancelled together with the infrared singularities at the level of the 
decay width, taking the bremsstrahlung process b — > s£ + £~ g into account. Gauge invariance 
implies that the QCD corrected matrix elements of the operators Oj can be written as 

{s£ + r\Oi\b) = F/ 9) (0 9 ) trce + PP(0 7 ) tiee , (5) 

where (0 9 ) tr ee and (C^tree are the tree-level matrix elements of 9 and 7 , respectively. 
Equivalently, we may write 

(Si + r \Oi\b) = [if ^O^ee + <0 7 ) tree ] , (6) 
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FIG. 1. Complete list of two-loop Feynman diagrams for b — ► 57* associated with the operators 
Oi and 02- The fermions (6, s and c quarks) are represented by solid lines, whereas the curly lines 
represent gluons. The circle-crosses denote the possible locations where the virtual photon (which 
then splits into a lepton pair) is emitted. 



where the operators O7 and 9 are denned as 



7 = ^0 7 ; 9 = ^0 9 . 

An 47T 



(7) 



We present the final results for the QCD corrected matrix elements in the form of eq. (§). 



A. Regularized 0{a s ) contribution of 0\ and 2 



The full set of the diagrams contributing to the matrix elements 

Mi={stT\Oi\b) (i = l,2) 



(8) 



at 0(a s ) is shown in Fig. ffl. As indicated in this figure, the diagrams associated with Ox 
and O2 are topologically identical. They differ only by the color structure. While the matrix 
elements of the operator O2 all involve the color structure 

AT2 _ 1 

y^ T a T a = c x c = _c 

a 

there are two possible color structures for the corresponding diagrams of Ox, viz. 

Tl = J2T a T b T a T b and r 2 = ^T a T b T b T a . 

a,b a,b 

The structure T\ appears in diagrams |l]a)-d) and r 2 in diagrams [Tje) and |I]f). Using the 
relation 
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I 



p 




J 



ap 



FIG. 2. The building blocks Ip and J a p which are used for the calculation of the two- loop 
diagrams |l|a)-d). The curly and wavy lines represent gluons and photons, respectively. 



we find that n = C T , 1 and r 2 



C T2 1 with 



N 2 c -l 
4iV c 2 



and C. 



4iV c 2 



Inserting N c = 3, the color factors are Cp 



|,C T1 = — | and C T2 = The contributions 
from 0\ are obtained by multiplying those from O2 by the appropriate factors, i.e. by 
C T1 /Cp = —1/6 and C T2 /Cp = 4/3, respectively. In the following descriptions of the 
individual diagrams we therefore restrict ourselves to those associated with the operator O2. 

In the current paper we use the MS renormalization scheme which is technically imple- 
mented by introducing the renormalization scale in the form /I 2 = /i 2 exp(7^) / (47r), followed 
by minimal subtraction. The precise definition of the evanescent operators, which is neces- 
sary to fully specify the renormalization scheme, will be given later. The remainder of this 
section is divided into 8 subsections. Subsections 1-6 deal with the diagrams |l]a)-d) which 
are calculated by means of Mellin-Barnes techniques ||451 . Subsection 7 is devoted to the 
diagrams |I|e) which are evaluated by using the heavy mass expansion procedure |46| . Among 
the diagrams |TJf ) only the one where the virtual photon is emitted from the charm quark line 
is non-zero. As it factorizes into two one-loop diagrams, its calculation is straightforward 
and does not require to be discussed in detail. It is, however, worth mentioning already at 
this point that it is convenient to omit this diagram in the discussion of the matrix elements 
of Oi and 2 and to take it into account together with the virtual corrections to 9 . Fi- 
nally, in subsection 8, we give the results for the dimensionally regularized matrix elements 
(sP-£-\Oi\b)(i = 1,2). 



1. The building blocks Ip and J a p 



For the calculation of diagrams pi) - d) it is advisable to evaluate the building blocks 
Jg and J a /3 first. The corresponding diagrams are depicted in Fig. After performing a 
straightforward Feynman parameterization followed by the integration over the loop mo- 
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mentum, the analytic expression for the building block Ig reads 



A 



■f^ r(c) /i 2e (1 - e) (r^ - r 2 7/3 ) ^ 3 X 



dx [x(l — x)] 



l-e 



r 2 - 



x(l — x) 



+ iS 



(9) 



where r is the momentum of the virtual gluon emitted from the c-quark loop. The term iS 
is the " ie-prescription" . In the full two-loop diagrams, the free index (3 will be contracted 
with the corresponding gluon propagator. Note, that Ip is gauge invariant in the sense that 
r% = 0. 

The building block J a p is somewhat more complicated. Using the notation introduced 



by Simma and Wyler ||47|| , it reads 



J, 



a;3 



16tt 2 



E(a, 0, r)Ai 5 + E(a, (3, q)Ai 6 - E(0, r, q) — Az 23 

q-r 

-E[a, r, g)^Az 25 - E(a, r, g)^Az 26 - E(0, r, q)^Ai 27 
q-r q-r q-r 



(10) 



where q and r denote the momenta of the (virtual) photon and gluon, respectively. The 
indices a and (3 will be contracted with the propagators of the photon and the gluon, 
respectively. The matrix E(a,f3,r) is defined as 



E(a, (3, r) = -{lalpf ~ flplcc) 
and the dimensionally regularized quantities Aik occurring in eq. fllOD read 

Ai 5 = AB + / dxdy [4(g-r) xy (1 — x)e + r 2 x (1 — x)(l — 2x)e 
Js 

+q 2 y(2-2y + 2xy-x)e + (1 - 3x)C] , 

Ai 6 = 4 5+ / cfe cfy [-4(g-r) ar y (1 - y)e - g 2 y (1 - y)(l - 2 y)e 

-r 2 x(2-2x + 2xy-?/)e- (l-3y)C] C~ 1_£ , 

-Ai 26 = 8 B + (q-r) dxdyxyeC~ 
Js 



(11) 



Ai 23 
Ai 25 



-l-e 



-8B (q-r) dxdyx(l — x)eC 
>s 



-l-e 



A* 27 = 8 5 + (g-r) /cfocfyy (1 - y) eC" 1 " 6 , 

where -B + = (1 + e)T(e) e 7Be /i 2e and C is given by 

C = m 2 — 2 x y(q-r) — r 2 x (1 — x) — q 2 y {I — y). 



(12) 
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The integration over the Feynman parameters x and y is restricted to the simplex S, i.e. 
y E [0,1 — x], x & [0, 1]. Due to Ward identities, the quantities Ai k are not independent of 
one another. Namely, 

q a J a /3 = and r 13 J a p = 
imply that Ai 5 and Ai s can be expressed as 

Ai 5 = Ai 23 H Ai 27 ; Ai e = — Ai 25 + Ai 26 . (13) 

q-r q-r 



2. General remarks 



After inserting the above expressions for the building blocks Jg and J a p into diagrams 
|l]a), b) and [L]c), d), respectively, and introducing additional Feynman parameters, we can 
easily perform the integration over the second loop momentum. The remaining Feynman 



parameter integrals are, however, non-trivial. In refs. |12] and f4~8 |, where the analogous 
corrections to the processes b — > S7 and b — > sg were studied, the strategy used to evaluate 
these integrals is the following: 

• The denominators are represented as complex Mellin-Barnes integrals (see below and 
refs. [Q,|U). 



• After interchanging the order of integration and appropriate variable transformations, 
the Feynman parameter integrals reduce to Euler f3- and V- functions. 

• Finally, by Cauchy's theorem the remaining complex integral over the Mellin variable 
can be written as a sum over residues taken at certain poles of (3- and V- functions. 
This leads in a natural way to an expansion in the small ratio z = m 2 /m 2 . 

However, this procedure cannot be applied directly in the present case: While the processes 
b — > 37 and b — >• sg are characterized by the two mass scales and m c , a third mass scale, 
viz. q 2 , the invariant mass squared of the lepton pair enters the process b — > s£ + £~. For 
values of q 2 satisfying 

q 2 q 2 

— 7? < 1 and < 1, 

mi Am 2 

most of the diagrams allow a naive Taylor series expansion in q 2 and the dependence of the 
charm quark mass can again be calculated by means of Mellin-Barnes representations. This 
method does not work, however, for the diagram in Fig. |T]a) where the photon is emitted 
from the internal s-quark line. Instead, we apply a Mellin-Barnes representation twice, as 
we discuss in detail in subsection 4. Using these methods, we get the results for diagrams 
|l]a)-d) as an expansion in s = q 2 jm\, z = m 2 /m 2 and s/(Az) as well as ln(s) and In (2). 
This implies that our results are meaningful only for small values of s. Fortunately, this is 
exactly the range of main theoretical and experimental interest in the phenomenology of the 
process b — > si + l~ . 



8 



3. Calculation of diagram ^b) 



We describe the basic steps of our calculation of the diagram in Fig. |l]b) where the 
photon is emitted from the internal 6-quark line. Our notations for the momenta are set up 
in Fig. |3|a). Inserting the building block Ip yields the following analytic expression for this 




FIG. 3. a) Momentum flow in diagram ||b), where the virtual photon is emitted from the 
6-quark line (see subsection 3); b) Momentum flow in the vertex correction diagram in Fig. |l|e) 
(see subsection 7). 



diagram: 



M 2 [lb] = %e ®*/ 8 C F r(e)e 27B y e (l - e)e i7re (47T] 



47T 2 

d d r 
(2?r) c 



dx 



[x(l — x)] 



l-e 



r 2 — m%/(x(l — x)) + i5] e 



v ' (p' + r) 2 — m£ (p + r) 2 — r 2 



Applying a Feynman parameterization according to 



1 



r(3 + 



D 1 D 2 D 3 Dl r(e) J s [ uDl + v d 2 + (i 



du dv dy y 



e-l 



u — v 



y)D 3 + yD,} 



3+e ' 



with 



D\ = (p' + r) 2 — m 



D 2 = (p + r) 2 - m 2 b , 
Di = r 2 — m 2 c /[x{\ — x)], 



D 3 = r\ 

and performing the integral over the loop momentum r, we obtain 

eQ d g 2 



M a [16] 



64tt 4 
i 



;i-e)C F r(2e)e 27B£ /i 4e x 



dx [x(l — x)} e / dvdudyy e u(p') 



Pi Pi 



(15) 



(16) 
(17) 



u{p), (18) 
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where the Feynman parameters u, v and y run over the simplex S, i.e u,v,y > and 
u + v + y < 1. P\, P2 and P3 are polynomials in the Feynman parameters and the quantity 
A;, reads 



Ab = ml(u + uv + v 2 ) — q 2 uv + 



ml y 
x(l — x) 



For q 2 < ml it is positive in the integration region. Therefore, one is allowed to do a naive 
Taylor series expansion of the integrand in q 2 . In order to simplify the resulting Feynman 
parameter integrals, it is convenient to first transform the integration variables x, y, u and 
v according to 



u 



(l-v')(-l+v' + u') 



(l-v')(l-u') 

; , x^x, y^yv. 

v' 



The integration region of the new variables is given by u' G [(1 — v'), 1] and v', x', y' G [0, 1]. 
Taking the corresponding Jacobian into account and omitting primes in order to simplify 
the notation, we find 



M 2 [lb] 



eQd 9 2 s 
64vr 4 



T -e)C F r(2e)e 27B y e x 



dx[x(l — x)\ e Jdv J du Jdy(vy) e (l-v)u(p') 

l-v 

where, in terms of the new variables, Ah reads 



Qi 



+ 



Q2 , Q 3 A fc 



u(p), (19) 



A b = ml(l - v)u + q' 



(l-v) 2 (l-v-u)(l-u) 



m„ 



vy 



'' x(l — x) 



Qi, Q2 and Q3 are rational functions in the new Feynman parameters. After performing the 
Taylor series expansion in q 2 , the remaining integrals are of the form 



dxdvdy / du[x(l — x)\ t (vy) 6 (1 — v) 



1 P(x,y,u,v) 



A 



l-v 



n+2e 
6,0 



(20) 



where P(x, y, u, v) is a polynomial in x, y, u, v; A^o = Ah(q 2 = 0); n and m are non-negative 
integers. We further follow the strategy used in [0,@ and represent the denominators A^ 
as Mellin-Barnes integrals. The Mellin-Barnes representation for [K 2 — M 2 )~ x reads (A > 0) 



(K 2 - M 2 ) x (K 2 ) x r(A) 2m 



ds 



M 2 

^2 



r(-s)r(A + s). 



(21) 



The integration path 7 runs parallel to the imaginary axis and intersects the real axis 
somewhere between —A and 0. The Mellin-Barnes representation of A^ is obtained by 
making the identifications 



K 2 <-> ml u{l — v ) and M 2 <-> — m 2 c y v/[x(l — x)]. 
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Interchanging the order of integration, it is now an easy task to perform the Feynman 
parameter integrals since the most complicated ones are of the form 

i 

J daa p{s \l - a) q ^ = (3(p(s) + 1, q(s) + l) . (22) 
o 

The integration path 7 has to be chosen in such a way that the Feynman parameter integrals 
exist for values of s 6 7. By inspection of the explicit expressions, one finds that this is 
the case if the path 7 is chosen such that Re(s) > — e. (Note that in this paper e is always 
a positive number). To perform the integration over the Mellin parameter s, we close the 
integration path in the right half-plane and use the residue theorem to identify the integral 
with the sum over the residues of the poles located at 

8 = 0, 1, 2, 3, 

s = 1 — e, 2 — e, 3 — e, 

s = l -2e, 2 - 2e, 3 - 2e, , 

8 = 1/2- 2e, 3/2 - 2e, 5/2 - 2e, (23) 

In view of the factor (m 2 c /m 2 ) s stemming from the Mellin-Barnes formula (|21|), the evaluation 
of the residues at the pole positions listed in eq. (^) corresponds directly to an expansion 
in z = m 2 c /m 2 . Note, that closing the integration contour in the right half-plane yields an 
overall minus sign due to the clockwise orientation of the integration path. After expanding 
in e, we get the form factors of M2[16] (see eq. (^)) as an expansion of the form 

F^\lb] = Y J ^l?z l \^{z), (24) 

i,l,m 

where i and m are non- negative integers and I is a natural multiple of | (see eq. (|2"3"D). 
Furthermore, the power m of ln(z) is bounded by four, independent of the values of i and 
/. This becomes clear if we consider the structure of the poles. There are three poles in 
s located near any natural number k, viz. at s = k, s = k — e and s = k — 2e. Taking 
the residue at one of them yields a term proportional to 1/e 2 from the other two poles. In 
addition, there can be an explicit 1/e 2 term from the integration over the two loop momenta. 
Therefore, the most singular term can be of order 1/e 4 and, after expanding in e, the highest 
possible power of ln(z) is four. 

4- Calculation of diagram [Jaj 

To calculate the diagram in Fig. [TJa) where the photon is emitted from the internal 
s-quark, we proceed in a similar way as in the previous subsection, i.e., we insert the build- 
ing block I j3, introduce three additional Feynman parameters and integrate over the loop 
momentum r. The characteristic denominator A a is of the form 

A a = (Am 2 b + Bq 2 + Cm 2 c +i6) 
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and occurs with powers 2e or 1 + 2e. The coefficients A, B and C are functions of the 
Feynman parameters. After suitable transformations, they read 



A = uv(l — v): B = uv z (l-u): C 



2„ .A. - 1/(1 



with it, x, y, v E [0, 1]. From this we conclude that the result of this diagram is not analytic 
in q 2 . We are therefore not allowed to Taylor expand the integrand. Instead, we apply the 
Mellin-Barnes representation twice and write 



1 l is I ^r( s + A)r(- s <)r (s <- s )e< 



AJ (BsTA Vy (2«) 2 r(A) 



Am 



fig 2 



Cmj 
Am 2 



(25) 



The integration paths 7 and 7' are again parallel to the imaginary axis and —A < Re(s) < 
Re(s') < 0. A takes one of the two values 2e and 1 + 2e. We have written eq. (|25| ) in such 
a way that non-integer powers appear only for positive numbers, i.e. we made use of the 
formula 

(x±t5) a = e ±ina {-xTiS) a . 

As in the previous subsection, the exact positions of the integration paths 7 and 7' are 
dictated by the condition that the Feynman parameter integrals exist for values of s and s' 
lying thereon. For A = 2e, we find that these integrals exist if 

-e < Re(s) < Re(s') < 0. 

Closing the integration contour for the s- and s'-integration in the left and right half-plane, 
respectively, and applying the residue theorem results in an expansion in s and z. As 
Re(s') > Re(s), the term T(s' — s) in eq. (|2~5"D does not generate any poles. For A = 2e, the 
poles which have to be taken into account are located at 

s' = 1 — e, 2 — e, 3 — e, s — — e, — 1 — e, — 2 — e, 

s' = 1 - 2e, 2 - 2e, 3 - 2e, s = -2e, - 1 - 2e, - 2 - 2e, 

s' = 0, 1, 2, 

For A = 1 + 2e, we find that the Feynman parameter integrals exist if 

-e < Re(s') < and - 1 - e < Re(s) < -2e . 

This condition implies that the poles at s — — e, — 2e in the above list must not be taken 
into account when applying the residue theorem. 

The final result for the form factors (eq.(^)) of this diagram is of the form 

if' 9) M = c 5,t^ lni (*V ln m W, (26) 

i,j,l,m 

where i, j, I and m all are non-negative integers. The remaining four diagrams in Fig. |l]a) 
and b) exhibit no further difficulties. 
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5. Calculation of diagrams [7|c ) 



Inserting the building block J a p allows us to calculate directly the sum of the two dia- 
grams shown in Fig. |I]c). After performing the second loop integral, one obtains 



M 2 [lc] = eQ n u £ C / (l + e) r(2e)e 2 ^y e e 2i - 



/dx dy du dv — ^— — - — rr— u(p' 



!l " " '" * 'u(p), (27) 



Pi P 2 P 3 A C 
+ — r + 



Al+2, A 2e A 2e 



where Pi, P 2 and P3 are polynomials in the Feynman parameters, which all run in the 
interval [0,1]. A c reads (using v' — v(l — u)) 

v' 

A c = m 2 b uv'y - q 2 yv' '(it + yv') - _ ^ [m 2 - q 2 y(l - x) (1 - y(l - x))\ . 

Note that we do not expand in q 2 at this stage of the calculation. Instead, we use the 
Mellin-Barnes representation (|2~TD with the identification 

K 2 <-> m 2 uv' y and M 2 ^ q 2 yv' (u + yv') -\ [m 2 — q 2 y(l — x) (1 — y(l — x))] . 

x[l — x) 

This representation does a good job, since (—M 2 /K 2 ) s turns out to be analytic in q 2 for 
s < Az, as in this range M 2 /K 2 is positive for all values of the Feynman parameters. We 
therefore do the Taylor expansion with respect to q 2 only at this level. Evaluating the 
Feynman parameter integrals as well as the Mellin-Barnes integral, we find the result as an 
expansion in z and s/(4z) which can be cast into the general form 



# 9) N = E c &l^^ln m W, (28) 



,,Lrn 



where i and m are non- negative integers and I = —i, — i + ~, + 1 



6. Calculation of diagrams ^d) 



After inserting the building block J a p and performing the second loop integral, the sum 
of the diagrams in Fig. [I]d) yields 



M 2 [ld] = 6 ^ff F (l + e) r(2e)e 2 ^V e 



256vr 4 



dx dy / dudv 



[x(l — x)] 



l+e 



u(p' 



Pi 



+ 



+ 



P 3 A d 



Af A 2 / j 



u(p), (29) 



where Pi, P 2 and P3 are polynomials in the Feynman parameters x, y, u and v. The 
parameters (x, y) and (u, v) run in their respective simplex. The quantity reads 



Ad — m h u I u 



yy 

1 — X 



+ q yv 



yv 



+ 



u _ (1 - y) 
1 — x x(l — x) 



m 2 v 
x(l — x) 
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Next, we use the Mellin-Barnes representation (^1|) with the identification 



K 2 <-> m\ u I u + 



1 — X 



q 2 yv 



yv 



u 



1-x) 



1 — x x(l — x) 



m 2 v 



x(l — x) 



Again, (—M 2 /K 2 ) s is analytic in q 2 for s < 4z, what allows us to perform a Taylor se- 
ries expansion with respect to q 2 . In order to perform the integrations over the Feynman 
parameters, we make suitable substitutions, e.g. 

(1- X >)(y> -(!-«')) 



x 



x , y 



i i 

U V , 



u 



u'(l-v'). 



(30) 



The new variables x', u', v' run in the interval [0, 1], while y' varies in [1 — v', 1]. Evaluating 
the integrals over the Feynman and Mellin parameters, we find the result as an expansion 
in z and sf(Az) which can be cast into the general form 



i,l,m 



(31) 



i and m are non-negative integers and I 



i + i, -2 + 1, 



7. Calculation of diagram \^e) 

We consider one of the diagrams in Fig. |l]e) in some detail and redraw it in Fig. |3|b). 
The matrix element is proportional to 1/A e , where 

A e = [(I - r) 2 - ml] [(l-q- r) 2 - m 2 ] [(/ - q) 2 - m 2 ] [I 2 - m 2 c ] r 2 . (32) 

q is the four-momentum of the off-shell photon, while / and r denote loop momenta. As 
q 2 < Am 2 in our application, we use the heavy mass expansion (HME) technique [[46] to 
evaluate this diagram. In the present CclSG, clS the gluon is massless, the HME boils down 
to a naive Taylor series expansion of the diagram (before loop integrations) in the four- 
momentum q. Expanding 1/A e in q, we obtain 

1 ST nt ' ' (q 2 ) l (q-ry(q-l) k 

- 2^ O e (n,m,i,j,k) 2 N(/ _ r)2 _ 2]m • 

n,m,i,j,k 1 cj lv ) cj 

Using the Feynman parameterization 

1 

1 T(n + m) f v m - % (l - v) n ~ l 



[I 2 - m 2 c ] n [(I - r) 2 - m 2 ] m T(n)T(m) J [I 2 - 2v(l ■ r) - m 2 + vr 2 ] n+m 



r v m - l {\-v) n - 1 

J [l 2 - 2v{l ■r)-m 2 + vr 2 } n + m ' ( ' J>i) 



o 



we can perform the integration over the loop momentum I. The integral over the loop 
momentum r can be done using the parameterization 



i 

i y r(i+p) f up- 1 



r 2 \ r 2 _ jnl \ ' Y(p) J ( 2 um g V+ 1 



du . (35) 
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The remaining integrals over the Feynman parameters u and v all have the form of eq. (|22|) 
and can be performed easily. The other two diagrams in Fig. |l]e) where the virtual photon 
is emitted from the charm quark can be evaluated in a similar way. The diagrams where 
the photon is radiated from the 6-quark or the s-quark vanish. 

As the results for the sum of all the diagrams in Fig. ^) are compact, we explicitly 
give their contribution to the form factors F a °' } (a = 1, 2; j = 7, 9). We obtain Fi 7) [le] = 0, 



ff[le] 



|F 2 (9) [le] and 



Ff[le] 



m c 



■lc 



124 
~27 



1 

e 

+ 



3 

12416 
3645 



128 
IE" 



s 

Tz 



s 



+ 



256 



105 

11072 
42525 



s 

Tz 



s 

~z 



2048 / 




+ 945 \ 




4971776 




4465125 


\Tz 



(36) 



8. Unrenormalized form factors of 0\ and O2 



We stress that the diagram |ljf ) where the virtual photon is emitted from the charm quark 
line is the only one in Fig. |l| which suffers from infrared and collinear singularities. As this 
diagram can easily be combined with diagram |]b) associated with the operator O9, we take 



it into account only in section [IV A| where the virtual corrections to O9 are discussed. 

The unrenormalized form factors Fa'^ of (s£ + £~\O a \b) (a = 1,2), corresponding to 
diagrams |l|a)-|l|e), are obtained in the form 



i,j,l,m 



m 'z) 



where i, j and m are non- negative integers and / = —i, — i + ~, — % + 1, We 

keep the terms with i and Z up to 3, after checking that higher order terms are small for 
0.05 < s < 0.25, the range considered in this paper. As we will give the full results for the 
counterterm contributions to the form factors in section prB| and the renormalized form 



factors in section [III Q and in appendix 0, it is not necessary to explicitly present the some- 
what lengthy expressions for the unrenormalized form factors. But, in order to demonstrate 
the cancellation of ultraviolet singularities in the next section, we list the divergent parts of 
the unrenormalized form factors: if> , if } , F 2 (7) and F 2 (9) : 
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F 2 9 L = ^ - orr^r ( 1890 + 1260 ™ + 5040 L » ~ 1260 L * + 252 S + 27 S 2 + 4 S 3 ) 



F (7) 
r 2Aiv 



+^T~ 420 + 2520 L M - 1260 L z + 2016 — + 1296 — + 1024 — 
2835 e \ V Az ) \ Az ) \Az 

92 



s \ „ f s \ . IS 



81 e' 

F idiv = -77^ + (1890 + 1260 m + 5040 - 1260 L s + 252 s + 27 S 2 + 4 S 3 ) 

' 243 76545 e v ' 

4 / , , „ / s \ „ I s \ is 



, -8085 + 2520 L, t - 1260 L z - 7056 — - 6480 — - 5_. , 
8505 e V " \Az J \Az J \Az 

(7) _ _ J6_ 

^l,div- 243e > 

where L s = ln(s), L 2 = ln(z), L M = In (^^j and z = m 2 /m 2 . 



B. C(a s ) counterterms to Oi and O2 

So far, we have calculated the two-loop matrix elements (s£ + £~\Ci Oi\b) (i = 1,2). As 
the operators mix under renormalization, there are additional contributions proportional to 
Cj. These counterterms arise from the matrix elements of the operators 

12 

SZijOj , 2 = 1,2, (38) 

3=1 

where the operators 0\-Ow are given in eq. (0). On and 12 are evanescent operators, i.e., 
operators which vanish in d = 4 dimensions. In principle, there is some freedom in the choice 
of the evanescent operators. However, as we want to combine our matrix elements with the 



Wilson coefficients calculated by Bobeth et al. [41], we must use the same definitions: 



On = (s Ll , lvla T a c L ) (c L ^YYT a b L ) -16 u (39) 
012 = {sLl»l v lvC L ) (cl^Yi^) - 16 2 . (40) 

The operator renormalization constants Zij = 5{j + 8Z{j are of the form 

O^o ( ni 1 11 \ f no 1 10 1 



Most of the coefficients a'™ needed for our calculation are given in ref. JIT]. As some are 



a 11 



not explicitly given in 


|41|) 


we list those for % = 


1, 2 and j = 


= 1,..., 


12: 




/-2 | -§ 


- 


16 n A 2 \ 12 _ 

27 U 12 9 \ "17 - 


58 12 
243 ' 19 


64 

729 ' 


„22 _ 

"19 ~~ 


1168 
243 


\ 6 00 | 000 


- 


•I 1 0/ ' 4 2 = 


116 _12 
81 ' u 29 ~~ 


776 

243 ' 


„22 _ 
u 29 ~~ 


148 
81 ' 



(42) 
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We denote the counterterm contributions to b — > s£ + £~ which are due to the mixing of 0\ 
or 2 into four-quark operators by i^Iq Uark and rk . They can be extracted from the 

equation 



loop 



Air 



,nct(7) 



->4quark 



(0 7 ) tree + Fl 



ct(9) 
+4quark 



<Og> 



tree 



(43) 



where j runs over the four-quark operators. As certain entries of a 11 are zero, only the 
one-loop matrix elements of 0%, 2 , O4, On and i2 are needed. In order to keep the 
presentation transparent, we relegate their explicit form to appendix 

The counterterms which are related to the mixing of Oi (i = 1, 2) into Og can be split 
into two classes: The first class consists of the one- loop mixing Oi — > Og, followed by taking 
the one- loop corrected matrix element of Og. It is obvious that this class contributes to 
the renormalization of diagram |l|f). As we decided to treat diagram |IJf) only in section 
IV A| (when discussing virtual corrections to Og), we proceed in the same way with the 
counterterm just mentioned. There is, however, a second class of counterterm contributions 
due to Oi — > Og mixing. These contributions are generated by two-loop mixing of 2 into Og 
as well as by one-loop mixing and one-loop renormalization of the g s factor in the definition 
of the operator Og. We denote the corresponding contribution to the counterterm form 
factors by F^P and F^g . We obtain 



F 



ct(9) 



22 12 

e 2 e 



a ig fig 



F 



ct(7) 







(44) 



where we made use of the renormalization constant Z„ s given by 



J 9* 



1 a s p 1 
4tt 2 e 



11 



-N f 
3 } 



N f = 5. 



(45) 



Besides the contribution from operator mixing, there are ordinary QCD counterterms. 
The renormalization of the charm quark mass is taken into account by replacing m c through 
Z mc ■ m c in the one-loop matrix elements of 0\ and 2 (see appendix |A]). We denote the 



corresponding contribution to the counterterm form factors by F i 
obtain 



ct(7) 



? ct(9) 
2,m c ren 



ct(7) _ pct(7) 
icven 2,m c ren 

32 



0; Fl 



1 105 + 84 — + 72 — 

945e \ \Az 1 1 ' 



ct(9) 



-F 



32 
2835 



105 + 1260 In -!— + — 336 + 1008 In ] + 



ct(9) 



2,m c ren 




S 

Tz 



396 + 864 In ^~ 

m c 



+ 



( 416 + 768 In — 

4z ) V m r 



(46) 
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where we have used the pole mass definition of m c which is characterized by the renormal- 
ization constant 

a s 4 /3 



Zm — 1 — 



4vr3 



(§ +( »»(£H 



(47) 



If one wishes to express the results for F^S 9 \ „ in terms of the MS definition of the charm 



quark mass, the expressions in eqs. ( f4"6| ) get changed according to 



where AF^ 9 l n „ reads 

t^llbcl Kill 



l,m c ren o 2,m c ren 



pct(9) _^ pct(9) »pct(9) 



,ct(9) 



(48) 



AE 



ct(9) 



2,m c ren 



64 
945 



105 + 84 



43 



+ 72 



s 

4z 



2 / - \ 3' 

S 



+ 64 



42 



, A* 2 
In — + - 

m r 3 



(49) 



We stress at this point that we always use the pole mass definition in the following, i.e., eqs. 

(E) for F'Zen- 

The total counterterms (i = 1, 2; j = 7, 9) which renormalize diagrams |l]a)-|l|e) are 

given by 

P ct (j) _ EiCtQ) pCt(j) pCt(j) /r n \ 

Explicitly, they read 



,ct(9) 



_ F ( 9 ) _ 



[5740 + 2520 rr 2 - 840 i 



2 2 - div 25515 

+ 840 L M (19 - 3 m - 54 L 2 + 48 L M ) + 3780 L 2 (-2 + 3 L 2 ) 

+ 420 L s (3 ztt + 2 + 6 L M ) - 630 L 2 + 252 3(1 - 2 L M ) - 54 L IM s 2 -2 s 3 (l + 4 L M ) 
+ 6048 (^) (18 L M - 9 L z - 1) + 7776 ( ^] (10 L M - 5 L z + 3) 



42 



42 



+ 1536 ( ^ ) (42L M -21L 2 + 19) 



F ct( 9 ) = _ F w^ + <_ [_ 62300 _ g40 m + 2520 n 2 



? (T) 

2, div 
,(9) 



2835 



(840 + 70 s + 7 s 2 + s 3 ) 



76545 



+ 840 L M (-3 irr - 54 L z + 48 L M - 791) + 3780 L 2 (3 L z + 88) 

+ 420 L s (3 w + 2 + 6 L M ) - 630 L 2 + s(252 - 504 L M ) - 54 s 2 L M - 2 s 3 (l + 4 L A 

- 6048 (-f) (28 + 90 L M - 45 L z ) - 7776 f — \ (27 + 62 — 31 L 2 ) 



42 



42 



- 768 



s 

Tz 



ct (7) _ p(7) 
^1 - _ ^l,div 



(295 + 564 — 282 L z 
1 



8505 



(840 L„ + 70 s + 7 s 2 + s 3 ) . 



(51) 
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The divergent parts of these counterterms are, up to a sign, identical to those of the un- 
renormalized matrix elements given in eq. (37), which proves the cancellation of ultraviolet 
singularities. 

As mentioned before, we will take diagram |I|f) into account only in section IV A. The 
same holds for the counterterms associated with the b- and s-quark wave function renor- 
malization and, as mentioned earlier in this subsection, the 0(a s ) correction to the matrix 
element of 5ZigOg. The sum of these contributions is 

a s a}} 



<^(Oi)l-loop + — [<^(0 9 )tree + (O 9 )i_ioop] , = Z^{m b )Z^{m s ) - 1, 

and provides the counterterm that renormalizes diagram [jj:). We use on-shell renormal- 
ization for the external b- and s-quark. In this scheme the field strength renormalization 
constants are given by 

^ W = 1 _f^(ii) 2< ('I + A +4 Y (52) 
47T 3 \mJ \e ei R / 

So far, we have discussed the counterterms which renormalize the 0(a s ) corrected matrix 
elements (s£ + £~\Oi\b) {i = 1,2). The corresponding one-loop matrix elements (of order 
O(a s )) are renormalized by adding the counterterms 

a s a i9 / n \ 

-. ^a/tree ■ 

47r e 



C. Renormalized form factors of 0\ and O2 



We decompose the renormalized matrix elements of Oi (i = 1, 2) as 



y(0) 



a 



,(9), 



(53) 



where 9 = ^ O g and O 



Oj. The form factors and F- 1 ' , expanded up to s 6 and 



(7) 



of the renormalized sum of diagrams |l]a)-e) read (L c = lnm c /mj = lnm c = L z /2) 



(9) 



1424 16 64 r , 

I in H L c ] L, 

729 243 27 



16 
243 



f. 



16 



32 



-1 



V1215 135 



+ 



2835 315 



L u s 2 + 



16 32 
76545 ~ 8505 



256 



- - ^ H ./:" • (54) 



-n ■ 



256 32 



+ 



128 

2 13 s] ~9 
16 



2 \ 77777 — 7J7 7T - L c ) + — L s + 



32 
51 



32 64 , . , 

+ — z- 1 )L U s 



405 45 



-2 



945 105 



32 



+ 



64 



25515 2835 
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The analytic results for f±, f[ 7 \ f!^\ and are rather lengthy. We decompose them as 
follows: 

ti b) = E <U 51 ^ zl L ™ + E A * L - ( 5? ) 

i,j,l,m i,j 

The quantities p^l? collect the half- integer powers of z = m 2 c /m1 = rh 2 c . This way, the 



summation indices in the above equation run over integers only. We list the coefficients 
K %im and Pah in appendix 0. 

If we give the charm quark mass dependence in numerical form, the formulas become 
simpler. For this purpose, we write the functions fa as 

fi b) = E k » H^') §l L * (a = 1, 2; 6 = 7, 9; 2 = 0, 3; j = 0, 1) . (58) 



i,3 



The numerical values for the quantities k a b \i,j) are given in Tab. [I] and [II] for m c = 0.25, 
0.29, 0.33. For numerical values corresponding to m c = 0.27, 0.29, 0.31 we refer to Tab. I 
and Tab. II in the letter version JO . 



Toy 



(9) 



„(9) 



(9) 



k\ 

k? 
k? 



I7J 



0.25 



-12.715 + 0.094699 i 
-0.078830 -0.074138 i 

-38.742 -0.67862 i 
-0.039301 - 0.00017258 % 

-103.83 - 2.5388 i 
-0.044702 + 0.0026283 i 

-313.75 - 8.4554 i 
-0.051133 + 0.022753 i 



0.29 



-11.973 + 0.16371 i 
-0.081271 - 0.059691 i 

-28.432 -0.25044 i 
-0.040243 + 0.016442 i 

-57.114 - 0.86486 i 
-0.035191 +0.027909 i 

-128.80 -2. 5243 i 
-0.017587 + 0.050639 i 



0.33 



-11.355 + 0.19217* 
-0.079426 -0.043950 i 

-21.648 - 0.063493 i 
-0.029733 + 0.031803 i 

-33.788 -0.24902 i 
-0.0020505 + 0.040170 i 

-59.105 - 0.72977 i 
0.052779 + 0.038212 % 



(7) 



,(7) 
"1 

,(7) 



k\ 



(7) 



.(7) 
"1 

.(7) 
"1 

,(7) 



-0.76730 - 0.11418 i 


-0.28480 - 0.18278 i 
-0.0032808 + 0.020827 i 
0.056108 -0.23357 i 
0.016370 + 0.020913 i 
0.62438 -0.027438 i 
0.030536 + 0.0091424 i 



-0.68192 - 0.074998 i 


-0.23935 -0.12289? 
0.0027424 + 0.019676 i 
-0.0018555 -0.17500 i 

0.022864 + 0.011456 i 

0.28248 - 0.12783 i 
0.029027 - 0.0082265 i 



-0.59736 - 0.044915 i 


-0.19850 - 0.081587 i 
0.0074152 + 0.016527 i 
-0.039209 - 0.12242 % 
0.022282 + 0.00062522 i 
0.085946 - 0.11020 i 
0.012166 -0.019772 i 



TABLE I. Coefficients in the decomposition of and for three values of m c . See eq. (|58|). 
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m c = 0.25 


m c = 0.29 


rh c = 0.33 


7 (9) 
2 


0. 


0) 


q 5D42 — 56819? 


fi RS38 — 98225 ?' 


4 S035 — 1 1 530 i 


i (9) 

«, 2 


;o. 


1) 


47298 4 D 44483 ?' 


48763 4- 35815 7 


47656 4 26370 ?' 


iA y ; 

2 


'i, 


o) 


7 4238 4- 4 0717 7 


3 3585 4 1 5026 i 


73780 4 38096 ?' 


7 (9) 
2 


'i, 


1) 


U.ZoOol T U.UUlUoOO ? 




u. i / o4u — u. lyuoz ? 


7 (9) 
2 


[2, 


0) 


33806 + 15 233 i 


-1 1906 + 5 1892 i 


-2.3570 4 1.4941? 


7 (9) 

LA ' 

2 


[2. 


1) 


26821 — 015770 i 


21115 — 16745 z 


012303 — 24102 ? 


7 (9) 
2 


'3. 


0) 


—42 085 4- 50 732 i 


—17 120 4 15 146 z 


—9 2008 4 4 3786? 


j (9) 

fi, 2 


;3. 


1) 


30680 — 1 3652 ?' 


10552 — 30383? 


—0 31667 — 22927? 

UiUlUU 1 V..' . — ■ — . .../ i_ 1 


, (7) 
2 




o) 


A fiD^S 4- f) 6851 i 


4 D91 5 -I- 44999 7 


q ^849 _i_ n 26949 i 




'0. 


1) 


n 

u 


n 


n 

VJ 


Jo 


1. 


o) 


1.7088 + 1.0967 i 


1.4361 4- 0.73732 i 


1.1910 4 0.48952? 


1.(7) 
2 


1. 


1) 


0.019685 - 0.12496 i 


-0.016454 -0.11806? 


-0.044491 -0.099160i 


too 

2 


% 


o) 


-0.33665 + 1.4014 i 


0.011133 4 1.0500?: 


0.23525 4 0.73452 i 


1.(7) 
2 


% 


1) 


-0.098219 -0.12548 i 


-0.13718 -0.0687337 


-0.13369 - 0.0037513 i 


1.(7) 
2 


% 


o) 


-3.7463 + 0.16463 % 


-1.6949 4 0.76698 i 


-0.51568 4 0.66118t: 


1.(7) 
«, 2 


3. 


1) 


-0.18321 - 0.054854 % 


-0.17416 4 0.0493597 


-0.072997 4 0.11863?: 



TABLE II. Coefficients in the decomposition of f 2 



(9) 



and /o 7 ^ for three values of rh c . See eq. 



IV. VIRTUAL CORRECTIONS TO THE MATRIX ELEMENTS OF 
THE OPERATORS 7 , O s , 9 AND O w 



A. Virtual corrections to the matrix element of Og and O 



'10 



As the hadronic parts of the operators Og and 0\q are identical, the QCD corrected 
matrix element of Oio can easily be obtained from the one of Og. We therefore present 
only the calculation for (s£ + £~\Og\b) in some detail. The virtual corrections to this matrix 
element consist of the vertex correction shown in Fig. ^b) and of the quark self-energy 
contributions. The result can be written as 



(s£ + £-\C 9 Og\b) = Cf ( ) F 9 w (0 9 ) tree + F^"(0 7 ) 



a 
4tt 



,(9), 



-(7), 



/tree 



(59) 



with Og = ?* Og and C, 



(o) 



In 



g(^9 



(0) 



a.ynr(l) 
4^°9 



We evaluate diagram §b) keeping the strange quark mass m s as a regulator of collinear 
singularities. The unrenormalized contributions of diagram |]b) to the form factors Fg 7 ^ and 
Fa read 
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IS 



<3 



_3_ 



b) 




d) 



-®- 



<2 

-3- 




FIG. 4. Some Feynman diagrams for 6 — > S7* or 6 — > si + l~ associated with the operators O7, 
Og and O9. The circle-crosses denote the possible locations where the virtual photon is emitted, 
while the crosses mark the possible locations for gluon bremsstrahlung. See text. 



2e / \ 2e 

tL\ ( _B 



r^(9) r.n \m b 1 4 \m b I 8 / A l^n 1 1 

F » I 4 " = " V-3 + ^T3 ( S + 2 S + 3 S + 2 

8 , , v 2 , 2 . , 16 20 A 16 . 2 116 

+ 3 (r)_ 3 (r) + y + y s + y ' 

F^[46] = -^(l + is + ^ 2 ) , (60) 

where we kept all terms up to s 3 . eiR and r = (ml /ml) regularize the infrared and collinear 
singularities in eq. (|60|). 

The b- and s-quark self-energy contributions are obtained by multiplying the tree level 
matrix element of O9 by the quark field renormalization factor 5Z^ = a/ Z^{va^)Z^{m^) — 1 , 
where the explicit form for Z^(m) (in the on-shell scheme) is given in eq. ([52j). 

Adding the self-energy contributions and the vertex correction, we get the ultraviolet 
finite results 

™ 16 20 A 16 ~ 116 „o , 

Fq y + y* + y* + yr* +/inf ' (61) 

/ \ 2e / \ 2e 

1 . , , ,|JL 



/inf = 7T (1 + s + - s 2 + - s 3 ) + ln(r) + - ln(r) - - ln 2 (r) , (62) 

eiR o zoo eiR o 6 

*$° = -h(i+b+?p) • (63) 



At this place, it is convenient to incorporate diagram |l|f ) together with its counterterms 



discussed in section [II B 
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It is easy to see that the two loops in diagram |TJf ) factorize into two one-loop contri- 
butions. The charm loop has the Lorentz structure of Og and can therefore be absorbed 
into a modified Wilson coefficient: The renormalized diagram ^f) is properly included by 
modifying in eq. (|59|) as follows: 



a 



(0) 



a 



(0,mod) 



^ 0) + 



a 



(0) 



+ 3 Cl 



(0) 



H n 



(64) 



where the charm-loop function H reads (in expanded form) 
1 



Hn 



2835 



I 200 2520 ln( — ) + 1008 I — ) + 432 ( — ) + 256 (— 



4z 



Az 



42 



(65) 



In the context of virtual corrections also the 0(e) -part of this loop function is needed. We 
neglect it here since it will drop out in combination with gluon bremsstrahlung. Note that 
H = h{z, s) + 8/9 ln(/x/m 6 ), with h defined in MM . 



B. Virtual corrections to the matrix element of O7 



We now turn to the virtual corrections to the matrix element of the operator 7 , con- 
sisting of the vertex- (see Fig. f|a)) and self-energy corrections. The ultraviolet singularities 
of the sum of these diagrams are cancelled when adding the counterterm amplitude 

a 7 

C 7 [Z 77 Z mb /Zl-l](s£ + £-\0 7 \b) trcc , where Z 77 = 1 - ^- . (66) 

The expressions for Z mb and Z gs are given in eqs. (|47|) and (flop , respectively. The renor- 
malized result for the contribution proportional to C 7 can be written as 

{stt \C 7 7 \b) = Cf (-^) [F 7 (9) (0 9 ) trce + F 7 (7) (0 7 ) tree ] , (67) 
with 7 = 7 and C 7 °^ = C 7 . The expanded form factors and read 

^ = -TH^r 2+ H' (68) 

„(7\ 32 32 r, 128 O . 

F} 7) = -L, + - + 8s + 6s 2 + — s 3 + f ini , (69) 

where the infrared- and collinear singular part f m { is identical to the one of Fg 9 ^ in eq. 
(p2f). Note that the on-shell value for the renormalization factor Z mb was used in eq. (|66|). 
Therefore, when using the expression for F 7 7,9 ^ in the form given above, the pole mass for 
rrib has to be used at lowest order. 
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C. Virtual corrections to the matrix element of 0$ 



Finally, we present our results for the corrections to the matrix elements of Os- The 
corresponding diagrams are shown in Fig. |]c) and d). Including the counterterm 

ot 16 

C 8 SZ 87 (s£ + £~\O r \b) tTce , where 5Z 87 = , 
yields the ultraviolet and infrared finite result 

(si + r \C 8 8 \b) = C® (-^) [F 8 (9) (0 9 ) trce + F 8 (7) (0 7 ) tree ] , (70) 

with = The expanded form factors and read 

r 9) 104 32 „ /1184 40 A A / 14212 32 „\ A « 

F 8 = — - — 7T 2 + — — - — 7T 2 )s + — — - — 7T 2 S 



9 27 



V 27 9 y V 135 3 J 

( 193444 560 2 \ 3 16 r , A „ 2 ^ 

+ 7T 2 s 3 H L s (1 + s + s 2 + s 3 ) , (71) 

V 945 27 J 9 V ; ' K ' 



(7) 32 8 2 44 8 A 2 40\ /32 2 316 ^ 2 

3 = L U H 71 ITT + — 7T S+ " 

s 9 M 27 9 9 V 3 3 7 V 9 9 



+ ^ i2 _^ j3 _8 i>(j + i2 + j3) , (72) 



V. BREMSSTRAHLUNG CORRECTIONS 



First of all, we remark that in the present paper only those bremsstrahlung diagrams 
are taken into account which are needed to cancel the infrared and collinear singularities 
appearing in the virtual corrections. All the other bremsstrahlung contributions (which are 
finite), will be given elsewhere [fP}| . 

It is known [p8 that the contribution to the inclusive decay width coming from the 
interference between the tree- level and the one- loop matrix elements of Og (Fig. |]b)) and 
from the corresponding bremsstrahlung corrections (Fig. f|f)) can be written in the form 



dr 



<)>) 



ds 

^99 

ds 



j-nvirt 
aL 99 

ds 



An 



+ 



^■pbrems 
Ul 99 



ds 



b,pole 



Kv tb \ 



48rr 3 



'l-SY(l + 2s) 



°9 



— UJg(s) 
7T 



(73) 



where C, 



(o) 



-Itt 

Os 



a 



(0) 



+ 4tt°9 



(1) 



This procedure corresponds to encapsulating the virtual 



and bremsstrahlung corrections in the tree-level calculation by replacing (Og) tree through 
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(l + — wg(s)) (Og)tree- The function ug(s) = o>(s), which contains all information on virtual 
and bremsstrahlung corrections, can be found in 



and is given by 



4 2 



2 2 5 + 4s 

-7T - 

9 



2s(l + s)(l-2s) 



3(1 -s) 



2s 



ln(s) + 



3(1 + 2s 
5 + 9s - 6s 2 
6(1 -s)(l + 2s) ' 



ln(l - s) 



(74) 



Replacing C® by Cf mod) (see eq. (@)) in eq. ([73]), diagram |]f) and the corresponding 
bremsstrahlung corrections are automatically included. 

For the combination of the interference terms between the tree-level and the one-loop 
matrix element of O-j (Fig. ^a)) and the corresponding bremsstrahlung corrections (Fig. 
) we make the ansatz 



dT 



77 



ds 
dT 77 
ds 



dT}f + dr^ cms 



ds 
4tt 



ds 

2 G 2 F m\ vole |TgT4h| 
48^3 



;i-s) 2 4(l + 2/s) 



a 



(0) 



0:.s 



7T 



U 7 (S) 



(75) 



where Cf^ = Cf~\ This time, the encapsulation of virtual and bremsstrahlung corrections 
is provided by the replacement {Octree — ► (l + ^^(s)) (C^tree- In order to simplify the 
calculation of uj^(s), we make the important observation that the form factors F^ and 
Fg 9 ^ have the same infrared divergent part f ini (eq. (BDp and (pT|)), whereas Fj 9 ^ and Fg 7 ^ 
are finite. Taking into account that in d dimensions the decay width dT(b — > s£ + £~)/ds 
corresponding to the matrix element 



M (b -> s£ + r) = (s£ + r |d} 0) of + or o 9 "> + oi» tree 



(o) n(o) 



(76) 



is given by 

dr(& -> Xg^+r) 

ds 

(l + (d-2)s) 



a; 



G 



(o) 



47T 
2 



2 Gi mi 



48tt 3 



;i-s) 2 (l + C(d-4)) x 



+ 



C 



+ 4(1 + (d-2)/s) 



a 



(0) 



+ 4(d 



y(o)p;(o)* 



l)Re[C^C£ 



, (77) 



one concludes that the combination 



Ar 



virt 



a 



(0) 



J-pvirt 
Ul 99 



(o) 



dY}f 



l + (d-2)s ds 4(l + (d-2)/s) ds 
is free of infrared and collinear singularities. Defining analogously 



(78) 



Ar 



brems 



a 



(0) 



jpbrems 
Ul gg 



(0) 



dr^ ems 



l + (d-2)s ds 4(l + (d-2)/s) ds 



(79) 
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and using the identity 



a 



(0) 



dTc 



a 



(0) 



dT 



77 



l + (d-2)s ds 4(l + (d-2)/s) ds 



Ar 



virt 



Ar 



brems 



(80) 



one concludes that also Ar brems is finite. This is because dTgg/ds and dT 77 /ds are finite due 
to the Kinoshita-Lee-Nauenberg theorem and because Ar virt is finite as mentioned above. 
The calculation of Ar brems is straightforward, as the integrand, expanded in e, leads to 
unproblematic integrals. Using the explicit results for Ar virt , Ar brcms and ujg(s), one can 
readily extract oo 7 (s) from eq. (|80|) : 



- In ' " 

3 \ m& 



U(s) 



9 



7T 



ln(s) ln(l - s) - 



1 8 + s 
3 2 + s 



ln(l 



2 1(2 - 2 s - s 2 ) , . , 1 16 -lis -17s 2 
m(s) — — 



3 1 



! (2 + S) 



18 (2 + s)(l-s) 



(81) 



The reasoning for the interference terms between the tree-level matrix element of 7 and the 
one-loop matrix element of Og and vice versa is analogous: We may combine this contribution 
with the corresponding bremsstrahlung terms coming from the interference of diagrams |]e) 
and f|f) making the ansatz 



dT 



79 



ds 

GTj79 

ds 



dT v 7g rt dTf 9 ems 



ds 



ds 



«em\ 2 GpTUf, le \ V t * s Vtb\ 



4rr 



48tt 3 



■(l-s) 2 12 



2Re (Cfcf* 



— uj 79 (s) 

71 



(82) 



The corresponding encapsulation is realized by the replacement (07,g)t 
(l + —u> 79 (§)) (Octree- This time, we make use of the fact that the quantities 



a pvirt 
mixed 



a 



(0) 



dT 



virt 



Re 



-1 



l + (d-2)s ds 
-2 



A(d-l) 



AT 



brems 
mixed 



a 



(0) 



brems 



Re 



'^(0)P?(0)* 
°7 °9 



f/s 



rfr brems 



l + (d-2)s ds 



4(d 



ds 



and 



(83) 
(84) 



are finite. For the function co»7 9 (s) we obtain 



^79 (sj 



-iln 
3 



2 s (3 - 2 s 



Li(s) 



7T 



9 (1-s) 



9 3 
1 5 -9s 

2 ln (s) + 7^T-ir- 



ln(s) ln(l - s) - 



1 2 + 7s 
9 s 



ln(l - s) 



(85) 



Note that the procedure described here does work only if one of the functions u) 7 (s), Ug(s) 
or w 7 g(s) is known already. 
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Finally, we remark that the combined virtual- and bremsstrahlung corrections to the 
operator 0\q (which has the same hadronic structure as 9 ) is described by the function 
u)g(s), too: 



dT 



10,10 



'/- 



jpvirt Jpbrems 
Ui 10,10 Ui 10,10 



dT 



10,10 



'/- 



ds 
ot em 

47T 



ds 



2 G%m 



b,pole \ VtsVtb 



48tt 3 



2s) 



(86) 



where C±q 



°io • 



VI. CORRECTIONS TO THE DECAY WIDTH FOR B -> X s £^£~ 



In this chapter we combine the virtual corrections calculated in chapters [TTT], [IV] and 
the bremsstrahlung contributions discussed in chapter [V| and study their influence on the 
decay width dT(b — > X s £ + £~)/ds. In the literature (see e.g. pi]), this decay width is usually 
written as 



rfr(6 -> X, 



Ait 



2 G 2 F m hpole \Vt s V t b\ 



(1 + 25) 



C, 



off 



+ 



48vT 3 
°10 



s) 2 X 



Cf + 12 Re (CfCf* ) ) , (87) 



+ 4(1 + 2/3) 

where the contributions calculated so far have been absorbed into the effective Wilson co- 
efficients Cj S , Cg S and . It turns out that also the new contributions calculated in the 
present paper can be absorbed into these coefficients. Following as closely as possible the 



'parameterization' given recently by Bobeth et al. f41 |) we write 



Ceff 
9 



cf 



°10 



7T 



u 9 (s) (A 9 + T 9 h(z, s) + U 9 h(l, s) + Wg h(0, s)) 



1 + 
1 + 



47T 

a s (jd) 



(c®f® + c®fW + a®fM 

a s (fi) 



(88) 



7T 



71 



-u,' 7 (.s) | A 7 - 

u 9 (s) ) A w 



4tt 



^(0)^(7) + C ^ IF ^ 



if i /.-') + ^ 0) F 8 (7) 



where the expressions for h(z,s) and u; 9 (s) (see eqs. ( p5|) and (|74D) were already available 
in the literature p4| , p8j^l| . The quantities ooi(s) and F^g, on the other hand, have been 
calculated in the present paper. We take the numerical values for A 7 , A 9 , A w , T 9 , Ug, and 
W 9 from fTf, while C{ 0) , and Af = cf ' eff) can be found in p8| . For completeness we 
list them in Tab. [TTT| . 

In Fig. |5| we illustrate the renormalization scale dependence of ReC^ (s). The dashed 
curves are obtained by neglecting the corrections calculated in this paper, i.e., ^(s), 
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u = 2 5 GeV 


u = 5 GeV 


u = 10 GeV 




0.267 


0.215 


0.180 




-0.697 


-0.487 


-0.326 


cf 


1.046 


1.024 


1.011 


V 7 ' 7 ) 


(-0.360, 0.031) 


(-0.321, 0.019) 


(-0.287, 0.008) 


Af 


-0.164 


-0.148 


-0.134 




(4.241, -0.170) 


(4.129, 0.013) 


(4.131, 0.155) 




(0.115, 0.278) 


(0.374, 0.251) 


(0.576, 0.231) 




(0.045, 0.023) 


(0.032, 0.016) 


(0.022, 0.011) 


«\ m 


(0.044, 0.016) 


(0.032, 0.012) 


(0.022, 0.009) 


(AO) Alh 


(-4.372, 0.135) 


(-4.372, 0.135) 


(-4.372, 0.135) 



TABLE III. Coefficients appearing in eq. (|88|) for [i = 2.5 GeV, [i = 5 GeV and \i = 10 GeV. 
For a s (/j,) (in the MS scheme) we used the two- loop expression with 5 flavors and a s (mz) = 0.119. 
The entries correspond to the pole top quark mass nit = 174 GeV. The superscript (0) refers to 
lowest order quantities while the superscript (1) denotes the correction terms of order a s . 



F 2 and F$ are put equal to zero in eq. flBlf ). The three curves correspond to the values 
of the renormalization scale /i = 2.5 GeV (lowest), /i = 5 GeV (middle) and fi = 10 GeV 
(uppermost). The solid curves are obtained by taking into account the new corrections. 
In this case, the lowest, middle and uppermost curve correspond to \i = 10 GeV, 5 GeV 
and 2.5 GeV, respectively. We conclude that the new corrections significantly reduce the 
renormalization scale dependence of ReCj G (s). 

Fig. |6] shows the renormalization scale dependence of ReCg ff (s). Again, the dashed 
curves are obtained by neglecting the new corrections in eq. (|88|) , i.e., and Fo 

are put to zero. We stress that wg(s) is retained, as this function has been known before. 
The three curves correspond to the values of the renormalization scale \i = 2.5 GeV (lowest), 
/i = 5 GeV (middle) and fi = 10 GeV (uppermost). The solid curves take the new corrections 
into account. Now, the lowest, middle and uppermost curve correspond to fi = 2.5 GeV, 5 
GeV and 10 GeV, respectively. We conclude that the new corrections significantly reduce 
the renormalization scale dependence of Re Cg S (s) , too. 

When calculating the decay width (0), we retain only terms linear in a s (and thus in uj, 
tug) in the expressions for (C^l 2 , |Cg ff | 2 and | Cicfl 2 - 111 ^ ne interference term Re (c^C^*^ 
too, we keep only linear contributions in a s . By construction, one has to make the replace- 
ments ujq —>■ ujjg and uj-j — > CJ79 in this term. 

Our results include all the relevant virtual corrections and those bremsstrahlung diagrams 
which generate infrared and collinear singularities. There exist additional bremsstrahlung 
terms coming e.g. from one-loop 0\ and O2 diagrams in which both, the virtual photon and 
the gluon, are emitted from the charm quark line. These contributions do not induce addi- 
tional renormalization scale dependence as they are ultraviolet finite. Using our experience 
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FIG. 5. The three solid curves illustrate the /j, dependence of ReCy ff (s) when the new correc- 
tions are included. The dashed curves are obtained when switching off these corrections. We set 
m c = 0.29. See text. 



from b — > 57 and b — > sg, these contributions are not expected to be large, but to give a 
definitive answer concerning their size, they have to be calculated . 



VII. NUMERICAL RESULTS FOR i? QU ARK(s) 



The decay width in eq. ( |37|) has a large uncertainty due to the factor m| pole . Following 
common practice, we consider the ratio 



R, 



quark (,"5 J 



1 dT(b -> X s £+t 



T(b^ X c ev) 



ds 



(89) 



in which the factor m\ le drops out. The explicit expression for the semi-leptonic decay 
width Y(b — > X c ev^ reads 



r(6 -> X c ev) 



m b,pole 
1927T 3 



V cb \ 2 g\ 



m 



c,pole 



K 



m 



b,pole 



i 1 I ' 

\ m b 



where g(z) = 1 — 8 z + 8 z 3 — z 4 — 12 z 2 \n(z) is the phase space factor, and 

= _ 2a s (m b ) fjz) 
V ' 3tt g(z) 



(90) 



(91) 
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FIG. 6. The three solid curves illustrate the \i dependence of ReC| ff (s) when the new correc- 
tions are included. The dashed curves are obtained when switching off these corrections. We set 
m c = 0.29. See text. 



incorporates the next-to-leading QCD correction to the semi-leptonic decay The func- 
tion f(z) has been given analytically in ref. |50f| : 



25 239 



25 
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z + — z 2 + z ln(z) 20 + 90 z - - z 2 + — z 



+. 2 ln 2 (.) (36 + z 2 ) + (1 - z 2 ) (11 - ^z + H hl (l - z) 

-4 (1 + 30 z 2 + z 4 ) ln(z) ln(l - z) - (1 + 16 z 2 + z 4 ) (6 U{z) - tt 2 ) 

1 - 



-32z 3/2 (l + 



7T" 



4 Li(v^) + 4 Li(— /£) - 2 hi(z) In 



1 + x/i 



(92) 



We now turn to the numerical results for -R qU ark(s) for 0.05 < s < 0.25. In Fig. |^ we 
investigate the dependence of -Rquark^) on the renormalization scale \i. The solid lines are 
obtained by including the new NNLL contributions, as explained in chapter [VI]. The three 
solid curves correspond to \i = 2.5 GeV (lowest line), \i = 5 GeV (middle line) and \i = 10 
GeV (uppermost line). The three dashed curves (again \i = 2.5 GeV for the lowest, fi = 5 
GeV for the middle and \x = 10 GeV for the uppermost line), on the other hand, show the 
results without the new NNLL corrections, i.e., they include the NLL results combined with 



the NNLL corrections to the matching conditions as obtained by Bobeth et al. 41 . From 
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FIG. 7. The three solid lines show the \i dependence of -Rquark(s) when including the corrections 
to the matrix elements calculated in this paper; the dashed lines are obtained when switching off 
these corrections. We set m c = 0.29. See text. 



this figure, we conclude that the renormalization scale dependence gets reduced by more 
than a factor of 2. Only for low values of s (s ~ 0.05), where the NLL /^-dependence is small 
already, the reduction factor is smaller. For the integrated quantity we obtain 

r0.25 

iWk = / dsR quaTk (s) = (1.25 ± 0.080*)) * 10~ 5 , (93) 

</0.05 

where the error is obtained by varying \i between 2.5 GeV and 10 GeV. Before our correc- 
tions, the result was -R qua rk = (1-36 ± 0.18) x 10~ 5 f4"Ifl . In other words, the renormalization 
scale dependence got reduced from ~ ±13% to ~ ±6.5%. 

Among the errors on -R qua rk('S) which are due to the uncertainties in the input parameters, 
the one induced by rh c = m c /mb is known to be the largest. We repeat at this point that m c 
enters (unlike in B — > X s ^) already the one- loop diagrams associated with 0\ and 0<i- We 
did the renormalization of the charm quark mass in such a way that m c has the meaning 
of the pole mass in the one-loop expressions. The meaning of m c in the corresponding 
two-loop matrix elements, on the other hand, is not fixed (for a discussion of this issue for 



B — > X s 7, see ref. [14]). As the running charm mass at a scale of 0{m^) is smaller than the 
pole mass, it numerically makes a difference whether one inserts a pole mass- or a running 
mass value for m c in the two-loop contributions. In a thorough phenomenological analysis 
this issue should certainly be included when estimating the theoretical error. We decide, 
however, to postpone the quantitative discussion of this point and will take it up when also 
the finite bremsstrahlung contibutions, which complete the NNLL calculation of -R quar k('S), 



are available [44]. For the time being, we interpret m c to be the pole mass in the two- loop 



contributions. In Fig. Rk) we vary rh c between 0.27 and 0.31, while in Fig. Rib) the more 
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FIG. 8. a) -Rquark(s) for m c = 0.27 (dashed line), m c = 0.29 (solid line) and rh c = 0.31 
(dash-dotted line) and fj,=5 GeV. b) i? qua rk('S) for rh c = 0.25 (dashed line), m c = 0.29 (solid line) 
and m c = 0.33 (dash-dotted line) and //=5 GeV. See text. 



conservative range 0.25 < rh c < 0.33 is considered. Comparing Fig. [7| with Figs. ||]a) and 
b), we find that at the NNLL level the uncertainty due to rh c is larger than the left-over 
//-dependence, even for the less conservative range of m c . For the integrated quantity -R qua rk 
we have an uncertainty of ±7.6% when rh c is varied between 0.27 and 0.31. Varying m c in 
the more conservative range, the corresponding uncertainty amounts to ±15%. 

A more detailed numerical analysis for -Rquarkf^) and -R qua rk, including the errors which 
are due to uncertainties in other input parameters as well as non-perturbative effects, will 
be given in ref. . 



To conclude: We have calculated virtual corrections of 0(a s ) to the matrix elements 
of Oi, 02, 07, 08, 09 and Oio- We also took into account those bremsstrahlung correc- 
tions which cancel the infrared and collinear singularities in the virtual corrections. The 
renormalization scale dependence of -R qua rk(s) gets reduced by more than a factor of 2. The 
calculation of the remaining bremsstrahlung contributions (which are expected to be rather 
small) and a more detailed numerical analysis are in progress |[44|| . 

Acknowledgements: C.G. would like to thank the members of the Yerevan Physics Insti- 
tute for the kind hospitality extended to him when this paper was finalized. We thank K. 
Bieri and P. Liniger for helpful discussions. 
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APPENDIX A: ONE-LOOP MATRIX ELEMENTS OF THE FOUR QUARK 

OPERATORS 



In order to fix the counterterms ^^J q u ark (i = 1,2) in eq. (f43"D , we need the one- loop 
matrix elements (s£ + t~\Oj\b)i.] oap of the four-quark operators Oi, 2 , 4 , On and Oi 2 . Due 
to the 1/e factor in eq. (0), they are needed up to O (e 1 ). The explicit results (in expanded 
form) read 



(s£ + r\0 2 \b) 1 . loop 



2c 



4 4 / / s \ / s 

— + -315 + 252 — + 108 

9e 2835 I V 4z 



2835 



I 05 rr" - 1 008 ' ^ ) + 128 ( ^ 



tree j 



64 



42 



(Al) 



(^ + r|o 1 |6) 1 _ loop = - (si + r|o 2 |6) Woop ^ 



(A2) 



( S £ + r|o 4 |6)i-ioo P 



m 6 



2c 



9 945 V ; 



(Or)t 



+ 



16 2 

+ (-420 + 1260 in - 1260 L s + 252 s + 27 s 2 + 4 s 3 ) 



27 e 8505 
4e 



+ (420ivr + 910 -630iL s vr -420L, -315tt 2 
8505 v 

+ 315 L 2 - 126 s + s 3 )] (Og) 



tree f , 



(A3) 



|On|6) 



l-loop 



64 / n 



27 \m, 



2c 



1 + 



4e / s 



+ 



5 \ 4z 7 V 4^ 



+ 



16 



63 V42 



(s£ + r|0 la |6) 1 . Ioop = - {st + V (Onl^Hoop 



<09> 



tree , 



(A4) 
(A5) 
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APPENDIX B: FULL s AND z DEPENDENCE OF THE FORM FACTORS Fj 



(7,9) 



In this appendix we give the dependence of fa (a = 1,2; b = 7,9) (see eq. (|57|)) on s 
and z. We decompose them as follows: 

fa'' = ^ K a}jlm ^ ^ ^ ^ ^ ^ ' 

i,j,l,m i,j 



The quantities p^ 6 ]- collect the half- integer powers of z = ml /ml = rh%. This way, the 
summation indices in the above equation run over integers only. On the following pages, we 
list the numerical values of Kaljim and for 

i — 0, 3; j = 0, 1; I — — 3, 3 and m = 0, 4. 

Coefficients not explicitly mentioned below vanish. 



Coefficients nfjj lm and p[ 9 ^ for the decomposition of 



(9) 
Pl,00 

(9) 

r 1,20 



3.8991 m d c 
-140.368 m r 



(9) 

o (9) 

PI, 30 



-23.3946 m c 

7.79821 m: 1 - 319.726 m r 



K 



(9) 

1,002m 



/ 



-4.61812 
14.4621 

-16.0864 
-14.73 



3.67166 i 
16.2155 i 
26.7517 i 
23.6892 i 



5.62963 
9.59321 
54.2439 
-28.5761 



I. 86168 i 

II. 1701 i 
14.8935 i 
34.7514 i 



-1.18519 
-15.4074 



7.44674 i 
29.787 i 








-0.790123 
-3.95062 



20.1481 



(9) 
K l,012m 



( 







-0.0493827 - 0.103427 i 
-0.592593 
4.95977 - 1.86168 i 
-9.20287 - 1.65483 i 



-1.18519 
-1.0535 



7.44674 i 
9.92898 i 









-2.37037 
3.16049 



K 



(9) 

l,10lm 



-2.48507 
4.47441 
71.3855 

-18.1301 
-72.89 



0.186168 i 
0.310281 i 
30.7987 i 
66.1439 i 
63.7828 i 



1.48148 
8.47677 
149.596 
-68.135 



1.86168 i 
33.5103 i 
67.0206i 
134.041 i 



12.5389 
-49.1852 



7.44674 i 
81.9141 i 








-0.790123 
-11.0617 



63.6049 








0.790123 





_(9) 

"1,112m 



/ 



-2.66667 
18.6539 
-41.6104 



1.86168 i 
7.44674 i 
3.72337 i 











-1.18519 

-4.74074 - 29.787 i -9.48148 

-2.37037 + 44.6804 i 14.2222 
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K 



(9) 

1,20/m 



/ 



-0.403158 
-0.0613169 
37.1282 
212.74 
-44.6829 
-137.203 





- 0.0199466 i 
+ 0.0620562 i 

- 1.36524 i 

- 52.2081 i 
+ 108.713 i 

- 106.832 i 









22.0621 - 1.86168 i 5.33333 0.790123 

-21.9215 - 52.1272 i 57.1724 - 7.44674 i -2.37037 2.37037 

272.015 - 163.828 i -119.111 - 156.382 i -21.3333 

-99.437 + 330.139 i 168.889 



K 0) 
^1,211™, 







0.0164609 
-5.33333 - 3.72337 i 
40.786 - 22.3402 i 
-111.356 











-2.37037 

-14.2222 - 67.0206 i -21.3333 

119.148 i 37.9259 



K 



(9) 



( 



-0.0759415 
0.00480894 
-1.81002 
79.7475 
425.579 
-87.8946 
-279.268 



0.00295505i 
0.00369382 i 
0.0871741 i 
1.72206 i 
76.6479 i 
148.481 i 
135.118 i 






-0.919459 

57.3171 - 1.86168? 
-68.8016 - 69.5029 i 

417.612 - 311.522 i 
-146.853 + 652.831 i 







-0.197531 

11.2593 2.37037 

129.357 - 7.44674 i -5.53086 4.74074 

-227.16 - 253.189 i -34.7654 

331.259 



K 



(9) 

l,3Hm 



/ 






















































0.0219479 
















-8.2963 - 5.58505 i 


-3.55556 











70.2698 - 49.6449 i 


-31.6049 


- 119.148 i 


-37.9259 








-231.893 + 18.6168 i 


11.8519 


+ 248.225 i 


79.0123 









\ 



Coefficients Kijj lm and p\'^ for the decomposition of j\ 



(7) 



(7) 
Pi, 00 

o {7) 

Pi, 20 



1.94955 nf c 
70.1839 m r 



o (7) 

Pi, 10 

o (7) 

Pl,30 



11.6973 m c 

-3.8991 C 1 + 159.863 m c 



(7) 



-1.14266 
-2.20356 
1.86366 
-1.21131 



0.517135 i 
1.59186 i 
- 3.06235 i 
+ 2.89595i 



+ 








-5.21743 + 1.86168 i 
-4.66347 



2.99588 



2.48225 i 











0.592593 + 3.72337 i 0.395062 

3.72337 i 0.395062 

-4.14815 



\ 



K 



(7) 

l,0Hm 



/ 



K 



(7) 

l,102ra 



-2.07503 + 1.39626i -0.444444 + 0.930842 i 



-25.9259 
11.4229 



+ 5.78065) 
- 15.2375) 



11.7509 + 15.6984) 



-3.40101 + 
-34.0806 + 
18.9564 - 



13.0318 i 
11.1701 i 
24.8225 i 











-4.4917 + 3.72337 i 0.395062 -0.395062 

10.3704 + 18.6168) 2.37037 

-14.6173 



\ 
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K 



(7) 

lAllm 



(7) 
K l,202m 



( 







-0.0164609 
1.03704 + 0.930842 i 

-4.66347 
,6.73754 + 1.86168 i 



( 






0.00555556 
-19.4691 + 1.59019i 
-90.4953 + 14.7788i 
23.8816 - 32.8021i 
38.1415 + 34.8683 i 












2 ^ 





















































0.592593 















7.44674 i 


2.37037 





) 






.18519 - 7.44674 i 


-2.37037 





o / 






















































-11.6779 + 0.930842 i 


2.96296 




-0.395062 






14.9329 + 22.3402i -24.438 + 3.72337i 1.18519 
-82.7915 + 39.0954 i 32.2963 + 44.6804 i 5.92593 
38.6436 - 80.673 i -41.5802 



-1.18519 





K 



(7) 

1,212m 



/ 







-0.0164609 
2.37037 + 1.86168 i 
-13.9904 + 3.72337 i 
27.5428 + 3.72337 i 











1.18519 

2.37037 + 22.3402 i 7.11111 

2.37037 - 29.787 i -9.48148 / 



K 



(7) 

l,30lm 



( 





0.00010778 + 0.00258567 i 
0.946811 - 0.0258567 i 
-41.9952 + 1.63673 i 
-189.354 + 25.8196i 
45.1784 - 52.4207i 
77.3602 + 54.2499i 






0.488889 

-30.2091 + 0.930842 i 

42.6566 + 31.0281 i 

-145.181 + 88.7403 i 

58.4491 - 184.927 i 







0.0987654 

-6.22222 -1.18519 

-57.765 + 3.72337i 2.76543 -2.37037 

70.9136 + 81.9141 i 11.0617 

-96.0988 



K 



(7) 

1,312m 
















































-0.0164609 














3.85185 + 2.79253 i 


1.77778 











27.3882 + 13.0318 i 


8.2963 + 44.6804 i 


14.2222 








69.4495 + 1.86168 i 


1.18519 - 74.4674 i 


-23.7037 









Coefficients ^Ijim an d P^Jj f° r the decomposition of 



{9) = -23.3946 m? 



P2,00 

o (9) 



= 842.206 rh r 



o {9) 

P2,10 

o {9) 

P2,30 



140.368 m c 

-46.7892 m; 1 + 1918.36 m c 



(9) 

2,002m 



K. 



(9) 

2,012m 



/ 














































































-24.2913 




22.0299 i 


-23.1111 


- 11.1701 i 















-86.7723 


+ 


97.2931 i 


-57.5593 


+ 67.0206 i 


7.11111 + 44.6804 


i 


4.74074 





96.5187 




160.51 i 


-325.463 


+ 89.3609 i 


92.4444 + 178.722 


i 


23.7037 





88.3801 


+ 


142.135i 


171.457 


- 208.509 i 


-120.889 




















































































0.296296 


+ 


0.620562 ? 



















3.55556 



















-29.7586 


+ 


11.1701 i 


7.11111 


+ 44.6804 i 


14.2222 











55.2172 


+ 


9.92898i 


6.32099 


- 59.5739i 


-18.963 





J 







\ 
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(9) 

2,102m 



0.8462 
-26.8464 
-428.313 
108.781 - 

437.34 + 



+ 1.11701 i 
+ 1.86168 i 
+ 184.792 i 
- 396.864 i 
382.697 i 







-8.88889 + 
-50.8606 + 
-897.575 + 402.124 i 
408.81 - 804.248 i 









11.1701 i 

201.062 i -75.2337 + 44.6804 i 4.74074 -4.74074 

295.111 + 491.485 i 66.3704 

-381.63 



K. 



(9) 

2,111m 



( 



16. 

-111.923 



11.1701 i 
44.6804 i 



249.663 + 22.3402 i 











7.11111 

28.4444 + 178.722i 56.8889 

14.2222 - 268.083i -85.3333 



\ 





K. 



(9) 

2,20lm 



( 





-0.0132191 + 0.11968 i 

0.367901 - 0.372337 i 

-222.769 + 8.19141 i -132.372 + 11.1701 i -32. -4.74074 

-1276.44 + 313.249 i 131.529 + 312.763 i -343.034 + 44.6804 i 14.2222 

268.098 - 652.279 i -1632.09 + 982.969 i 714.667 + 938.289 i 128. 

823.218 + 640.989 i 596.622 - 1980.83 i -1013.33 








-14.2222 





(9) 
K 2,2llm 



( 







-0.0987654 
32. + 22.3402 i 
-244.716 + 134.041i 
668.137 











14.2222 

85.3333 + 402.124 i 128. 

-714.887 i -227.556 



K. 



(9) 

2,302m 



0.0142243 + 0.0177303 i 

0.0288536 - 0.0221629 i 

10.8601 - 0.523045 i 

-478.485 + 10.3323 i 

-2553.47 + 459.887 i 

527.368 - 890.889 i 

1675.61 + 810.709 i 






5.51675 
-343.902 + 11.1701 i 

412.809 + 417.017 i 
-2505.67 + 1869.13 i 

881.117 - 3916.98 i 






1.18519 
-67.5556 
-776.143 + 44.6804) 
1362.96 + 1519.13) 
-1987.56 







-14.2222 
33.1852 
208.593 









-28.4444 





K. 



(9) 

2,312m 



/ 
















































-0.131687 














49.7778 + 33.5103 i 


21.3333 











421.619 + 297.87 i 


189.63 + 714.887) 


227.556 








1391.36 - 111.701 i 


-71.1111 - 1489.35) 


-474.074 









Coefficients i^ljim an d Pzjj f° r the- decomposition of 

pW Q = - 11.6973 ml pgo = - 70.1839 m c 

P2T20 = - 421.103 m c pgo =23.3946 m' 1 - 959. 179 m c 
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/ 



6.85597 
13.2214 
-11.182 
7.26787 







+ 3.10281i 

- 9.55118 i 
+ 18.3741) 

- 17.3757 i 








31.3046 - 11.1701) 
27.9808 
-17.9753 + 14.8935) 








-3.55556 - 22.3402 i 
-22.3402 i 
24.8889 








-2.37037 
-2.37037 




/ 



12.4502 
155.555 
. -68.5374 
V-70.5057 



8.37758i 
34.6839 i 
91.4251 i 
94.1903) 



2.66667 
20.4061 
204.484 
-113.738 



5.58505) 
78.1908i 
67.0206* 
148.935) 



26.9502 
-62.2222 



22.3402 i 
111.701* 



87.7037 








-2.37037 
-14.2222 









2.37037 









0.0987654 
-6.22222 - 5.58505 i 

27.9808 
-40.4253 - 11.1701 i 



\ 







-3.55556 

-44.6804 i -14.2222 

-7.11111 + 44.6804) 14.2222 



/ 



-0.0333333 

116.815 - 9.54113 i 70.0677 - 5.58505) 

542.972 - 88.6728) -89.5971 - 134.041) 

-143.29 + 196.813) 496.749 - 234.572) 

V-228.849 - 209.21) -231.862 + 484.038) 









17.7778 2.37037 

146.628 - 22.3402) -7.11111 7.11111 

-193.778 - 268.083) -35.5556 

249.481 



\ 







0.0987654 
-14.2222 - 11.1701 i 
83.9424 - 22.3402) 
-165.257 - 22.3402) 











-7.11111 

-14.2222 - 134.041) -42.6667 

-14.2222 + 178.722) 56.8889 



\ 



0.000646678 - 0.015514) 
-5.68087 + 0.15514) 
251.971 - 9.82039) 
1136.13 - 154.918) 
-271.07 + 314.524) 
-464.161 - 325.499) 






-2.93333 
181.255 - 5.58505) 
-255.94 - 186.168) 
871.089 - 532.442) 
-350.695 + 1109.56) 







-0.592593 

37.3333 7.11111 

346.59 - 22.3402) -16.5926 14.2222 

-425.481 - 491.485) -66.3704 

576.593 



/ 



V- 






















































0.0987654 
















23.1111 - 16.7552) 


-10.6667 











164.329 - 78.1908 i 


-49.7778 


- 268.083) 


-85.3333 








416.697 - 11.1701 i 


-7.11111 


+ 446.804) 


142.222 
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